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Abstract 

The strong coupling constant is one of the fundamental parame¬ 
ters of the standard model of particle physics. In this review I will 
briefly summarise the theoretical framework, within which the strong 
coupling constant is defined and how it is connected to measurable 
observables. Then I will give an historical overview of its experimen¬ 
tal determinations and discuss the current status and world average 
value. Among the many different techniques used to determine this 
coupling constant in the context of quantum chromodynamics, I will 
focus in particular on a number of measurements carried out at the 
Large Electron Positron Collider (LEP) and the Large Hadron Col¬ 
lider (LHC) at CERN. 
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1 Introduction 


The strong coupling constant, a s , is the only free parameter of the lagrangian 
of quantum chromodynamics (QCD), the theory of strong interactions, if 
we consider the quark masses as fixed. As such, this coupling constant, or 
equivalently g s = \[Afiux s , is one of the three fundamental coupling constants 
of the standard model (SM) of particle physics. It is related to the SU(3)c 
colour part of the overall SU(3)c x SU(2)^ x U(l)y gauge symmetry of the 
SM. The other two constants g and g' indicate the coupling strengths relevant 
for weak isospin and weak hyper charge, and can be rewritten in terms of 
the Weinberg mixing angle tan$w = g'/g and the fine-structure constant 
a = e 2 /(47r), where the electric charge is given by e = gsmOyj. Note that 
natural units are used throughout. 

While typically denoted as constants, actually all these coupling strengths 
vary as a function of the energy scale or momentum transfer Q of the par¬ 
ticular process looked at, as will be discussed later. In contrast to a(Q 2 ), 
which increases with increasing Q, the strong coupling a s (Q 2 ) decreases for 
increasing scale, leading to the famous property of QCD known as asymp¬ 
totic freedom. It is interesting to compare the values of these two coupling 
strengths at some fixed scale, such as the mass of the Z boson, Q cs M z ~ 91 
GeV. We find that a/M|) ~ 1/128, whereas ay(M|) ~ 0.12; that is, the 
strong coupling is still about 15 times larger than the fine-structure constant 
at energy scales much larger than those relevant for quark confinement into 
hadrons (Q ~ 1 GeV). Thus, strong interactions are indeed strong com¬ 
pared to electroweak interactions, even at large energy scales such as those 
probed by CERN’s past and present colliders, in particular the Large Elec¬ 
tron Positron Collider (LEP) or the Large Hadron Collider (LHC). 

The different energy dependence of the coupling strengths triggers the im¬ 
mediate question if and at which exact energy scale these coupling constants 
become of equal strength, implying the onset of a possible grand unification. 
Obviously, the answer to this question also depends on the precision at which 
a and a s have been determined by experiment, and it is instructive to realize 
that today a(Q 2 = 0) ~ 1/137 is known at an accuracy of 32 parts per bil¬ 
lion [T], whereas the relative uncertainty of the current world average (WA) 
value [2] of a s (M|) = 0.1185 ± 0.0006 amounts to half a percent; quite an 
astonishing difference. 

Besides the wish to improve the accuracy of the aforementioned very 
high energy extrapolation, it is of general importance to know a s at the best 
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possible precision, since it enters the calculation of each and every process 
that involves strong interactions and thus ultimately limits the precision 
at which such processes can be predicted theoretically. As a most recent 
and prominent example, it is worth mentioning that the uncertainty on a s 
gives a non-negligible contribution to the overall theoretical uncertainty on 
Higgs boson production at the LHC [3]. Correspondingly, this limits the 
studies looking for effects beyond the SM that could manifest themselves 
through deviations of the measured Higgs production cross sections from 
their theoretical predictions. In the following I will indicate the experimental 
and theoretical difficulties that limit the precision at which we know this 
fundamental parameter, but also highlight the dramatic improvements, which 
have been achieved during the last three decades. 

2 Theoretical framework 

The basic elements of QCD, including a discussion of the scale dependence 
of a s (Q 2 ) and the related structure of theoretical predictions obtained in 
perturbation theory, are summarized elsewhere in this series of reviews [TJ. 
Further extensive descriptions of the theoretical framework can be found 
in, e.g., Refs. mun. Here I will only highlight a few important aspects 
of perturbative QCD (pQCD), that are relevant for the remainder of this 
review. 

When calculating amplitudes corresponding to Feynman graphs that in¬ 
volve loop diagrams, ultraviolet divergences are encountered. The procedure 
of renormalization absorbs these divergences into a redefinition of the bare 
parameters and fields that appear in the lagrangian. In particular, this leads 
to the renormalised or so-called running coupling constant a s (/i 2 ), a function 
of the (unphysical) renormalization scale /i. If /i is chosen close to the scale of 
the momentum transfer Q in a given process, then a s (/i 2 ~ Q 2 ) is indicative 
of the effective strength of the strong interaction in that process [2j. This ex¬ 
plains why in the literature we often find a discussion of the running coupling 
constant as function of the physical scale Q , while the renormalized coupling 
actually is a function of the unphysical scale fi. This will also become clearer 
from the following discussion of the structure of perturbative predictions. 

While the value of a s (/i 2 ) at a fixed scale /i can not be predicted and has 
to be determined from experiment instead, its fi dependence is given by the 
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renormalization group equation, 


= PM = ~ b o a2 s + b i a l + ®( a t) 


( 1 ) 


The first two coefficients in the perturbative expansion of the so-called /3- 
function of QCD are b 0 = (33 — 2n/)/(127r) and b\ = (153 — 19nj)/(247r 2 ), 
where rif is the number of “light” quark flavours (m q <C /i). Most impor¬ 
tantly, for rif < 17 we have b 0 > 0, which leads to a decreasing coupling 
strength for increasing scale (asymptotic freedom), as originally predicted 
by Politzer [7], Gross and Wilczek [8j. Considering only the first term of 
the expansion on the right hand side of eq. [lj a solution can be written as 
« s (h 2 ) — (6ol n (A* 2 /A 2 )) - \ with A ~ 200 MeV defined as the scale where 
a s (/i 2 ) formally diverges. Whereas at small scales of order GeV or lower the 
coupling constant increases dramatically, making any perturbative approach 
to the solution of low-energy strong interactions and the property of con¬ 
finement meaningless, it is the property of asymptotic freedom that leads to 
an expansion parameter a s well below unity and thus allows perturbative 
methods to be applied for the calculation of scattering processes. 

To second order, including the resummation of leading logarithms of type 
In (/j, 2 /Q 2 ), a solution of eq. [l] allows to relate a s at one scale /r 2 to that at 
another scale Q 2 , 



In the following we will use the resulting expansion 



(3) 


This shows that a change of scale only manifests itself as a non-leading effect 
in a s ; in other words, a meaningful determination of the running coupling 
constant necessarily has to involve a next-to-leading order (NLO) prediction. 
In order to highlight this even further, let’s look at the perturbative structure, 
up to NLO, of some generic cross section that is proportional to a s at leading 
order (e.g., a three-jet cross section in e + e _ annihilations at = Q), 



+ 0(a 3 s ) . (4) 
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The coefficients A and B have to be calculated for the specific process at 
hand. Now let us first assume that only A is known for a particular process 
”1”, i.e., only the leading order (LO) expansion is available, a\° = a s (/i 2 )Hi. 
However, at the same LO we could equally well write this prediction as 
a\° = a s ([i 2 )Ai = a s (Q 2 )A 1 , because using the above expansion of the cou¬ 
pling constant, eq. [3j we see that the scale dependence would only appear as 
an NLO correction, namely a\°' = a s (Q 2 )Ai — a 2 (Q 2 )boAi \n/i 2 /Q 2 . Thus, 
strictly sticking to the LO expression, it is clear that an experimental mea¬ 
surement of o\ and its comparison to cr]" 0 only allows to determine some 
“effective” LO coupling constant ct® 5,1 , where it is unclear to which scale 
this really corresponds to. Furthermore, repeating the same procedure at 
LO for some other process ”2”, at some different physical energy or momen¬ 
tum scale, would result in a measurement and most likely these two 

measurements of the effective coupling constant would give differing results. 

Looking again at the expression a\° = a s (Q 2 )Ai — a 2 (Q 2 )b 0 A 1 In p 2 /Q 2 
we also see that a\° strongly depends on the unphysical scale p, since the 
logarithm with the explicit p dependence already appears at NLO. Corre¬ 
spondingly, a determination of a s (Q 2 ) using this prediction would result in 
a large uncertainty when varying the unphysical parameter p in the fits to 
the measured cross section. This procedure of p-variations, typically over 
a range of 0.5 < n/Q < 2, has become a standard approach to estimat¬ 
ing the possible impact of unknown higher-order contributions. In fact, the 
p-dependence always appears at one order higher than the fully known per¬ 
turbative expansion. More concretely, let’s now assume that also the NLO 
coefficient B has been calculated. Then, by plugging the expansion [3] into 
eq. [4] we find 



We see that there is no p-dependence up to NLO; at this order it is thus 
equivalent to set fi = Q and to write a(a s (fi 2 ), ^) = o‘(«s(Q 2 )); be., we can 
replace the dependence of the running coupling constant on the unphysical 
scale p with a dependence on the physical scale Q of the process at hand. 
Furthermore, we see that the explicit p-dependence of the cross section pre¬ 
diction only appears at next-to-NLO (NNLO), i.e. suppressed by two powers 
of a s relative to the LO term. This leads to a smaller uncertainty of the 
extracted a s (Q 2 ) value when varying p in the fit procedure. Finally, the 
NLO expression in eq. [5] leads to the first non-trivial dependence of the cross 
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section on a s (Q 2 ) at the particular scale Q. Therefore, determinations of 
a s (Q 2 ) from two different processes, at possibly different values of Q , using 
the NLO predictions for the cross sections and the running of a s in order to 
relate the measured values to each other, should result, within uncertainties, 
in consistent measurements. Similarly, the value of a s (Q 2 ) determined at 
NLO for some specific process can be used for predicting, at the same order, 
the cross section for another process at a different physical scale. 

The extension of this discussion to NNLO and beyond is straightforward, 
and easily shows that theoretical uncertainties estimated from p-variations 
should decrease even further. This is nicely illustrated in Fig. [lj where the 
dependence of the extracted value of a s (M|) is shown, when using the LO, 
NLO and NNLO pQCD expressions for fitting the measured hadronic width 
of the Z boson, normalised to its leptonic width [9j. Ultimately, for an 
observable known at all orders in pQCD the /i-dependence would vanish 
completely, as it should be for a physical observable that cannot depend 
on unphysical parameters. In fact, the p-dependence of the NLO term in 
expression [4] could have simply been derived by imposing this requirement 
for a physical observable and using the renormalization group equation. 

At this stage it should have become clear, but still is worth highlighting, 
that the running coupling constant a s (Q 2 ) is not a physical quantity, but sim¬ 
ply a parameter defined in the context of a particular theoretical framework, 
namely pQCD up to some order in ot s . It can be determined experimentally 
in this context and used for making predictions for observables calculated 
within the same framework. 

I would like to conclude these theoretical considerations by highlighting a 
further consequence of the particular scale behaviour of a s : An uncertainty 
5 on a measurement of a s (Q 2 ), at a scale Q, translates to an uncertainty 
S' = (a 2 (M%) / a 2 (Q 2 )) ■ 5 on a s (M|); that is, S' < S for Q < M z . This 
enhances the impact of precise low-Q measurements, such as from r decays 
(c.f. below), in combinations performed at the M z scale. 

3 Observables 

The strong coupling constant has been measured in a large variety of physics 
processes, using many different observables. As depicted in Fig. [2| sensitiv¬ 
ity to the coupling of quarks to gluons is obtained by studying, e.g., deep- 
inelastic lepton-nucleon scattering, e + e _ annihilations, hadron collisions or 
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Figure 1: ot s (M |) determined from the scaled hadronic width of the Z 
boson, in LO, NLO and NNLO QCD, as a function of the renormalization 
scale = yu/M z ; taken from Ref. pi]. 

resonance decays. Since we are not able to directly measure partons (quarks 
or gluons), but only hadrons and their decay products, a central issue is estab¬ 
lishing a correspondence between observables obtained at the partonic and 
the hadronic level. The only theoretically sound correspondence is achieved 
by means of infrared and collinear safe (ICS) quantities (see e.g. Ref. [2]), 
which allow to obtain finite predictions at any order of perturbative QCD. 
ICS observables are insensitive to the addition of a soft parton, or to the 
splitting of one parton into two collinear ones. This guarantees that singu¬ 
larities, which appear in the infrared or collinear limits of diagrams involving 
real and/or virtual radiation, cancel order by order in perturbation theory. 

Generally speaking, ICS observables can be divided into different classes, 
depending on how detailed the hadronic final state is being scrutinized. As for 
example discussed in Ref. [2], the simplest case of ICS observables are total 
cross sections. More generally, when measuring fully inclusive observables, 
the final state is not analyzed at all regarding its (topological, kinematical) 
structure or its composition. Basically the relevant information consists in 
the rate of a process ending up in a partonic or hadronic final state. 
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Figure 2: Examples of Feynman diagrams describing hadronic final states 
in processes which are used to measure a s ; taken from Ref. [9]. 


In e + e~ annihilation, widely used examples are the ratios of partial widths 
or branching ratios for the electroweak decay of particles into hadrons or lep¬ 
tons, such as Z or r decays. The strong suppression of non-perturbative 
effects, (P(A 4 /Q 4 ), is one of the attractive features of such observables. How¬ 
ever, at the same time the sensitivity to radiative QCD corrections is small, 
since here the perturbative expansion is of the type 1 + c n a 7> correspond¬ 
ing to, e.g., a 4% correction, 1+<x,(M|)/7T ~ 1 + 0.04, for the scaled hadronic 
Z width. In the case of r decays not only the hadronic branching ratio is 
of interest, but also moments of the spectral functions of hadronic tau de¬ 
cays, which sample different parts of the decay spectrum and thus provide 
additional information. 

Other examples of fully inclusive observables, that are used for a s determi¬ 
nations, are the total top-pair production cross section in proton-proton colli¬ 
sions, the ratio of the hadronic to leptonic branching ratios of T decays, which 
is proportional to at LO (cf. the right-most diagram in Fig. [2|, or struc¬ 
ture functions (and related sum rules) in deep-inelastic scattering. Struc¬ 
ture functions are sensitive to a s through the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi unnnniM evolution equations, e.g. dF 2 (x, Q 2 )/dlnQ 2 oc 
a s (Q 2 ) P qg <g) g(x, Q 2 ), which depicts in a simplified manner the scaling viola¬ 
tion of the F ‘2 structure function; here x is the Bjorken scaling variable, P qg is 
a so-called splitting function and g(x, Q 2 ) is the parton distribution function 
(PDF) of the gluon. Such equations are used in global PDF fits in order to 
relate measurements at different Q scales to each other and to fit the PDFs 
at a chosen initial scale. We see that in such approaches the fit results for 
a s and g(x, Q 2 ) are strongly correlated. Similar considerations apply to the 
measurements of scaling violations of fragmentation functions. 

Compared to inclusive observables, the particular structure, topology or 
composition of the hadronic final state can give enhanced sensitivity to a s , 
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therefore cross sections differential in one or more variables characterizing 
this structure are of interest. The obvious example is the measurement of jet 
cross sections and jet rates, where the probability of producing an additional 
jet is directly proportional to a s (for a general discussion of jets and jet 
algorithms I refer the reader to, e.g., Refs. 013 and references therein). 
Besides jet quantities, another class of observables, so-called event shapes , 
have been widely used, in particular for measurements in e + e~ annihilations, 
but also in lepton-proton and hadron collisions. The classic example of an 
event shape is the Thrust [Q, [15] (T or f = 1 — T) in e + e _ annihilations, 
defined as 


T 


max 

n T 


E» I Pi • n r 
Ei \Pi\ 


( 6 ) 


where /7, are the momenta of the particles or the jets in the final-state and 
the maximum is obtained for the Thrust axis n T . In the Born limit of the 
production of a perfect back-to-back quark-antiquark pair the limit f —> 0 is 
obtained, whereas a perfectly symmetric many-particle configuration leads to 
r —> 1/2. Figure [3] (left) shows an example of measurements by the ALEPH 
experiment at different centre-of-mass energies. 

Besides Thrust, many other similar observables such as C-parameter, 
Heavy Jet mass, Jet Broadening or the differential three-jet rate were pro¬ 
posed and used for a s -determinations. They all provide a measure of the 
topology of an event, and typically are defined such that they take on small 
values for pencil-like (back-to-back) configurations, and large values for more 
spherical topologies that arise from single or multiple hard gluon radiation. 
This provides sensitivity to a s at LO in perturbation theory, with normal¬ 
ized cross sections expressed as an expansion of the type eq. [4] As discussed 
further below, predictions are known up to NNLO and complemented by the 
all-orders resummation of logarithms of the event-shape variable (i.e., terms 
of the form a n s In'" r). An important aspect of event-shape variables is their 
enhanced sensitivity to non-perturbative effects compared to more inclusive 
quantities, with power corrections of ~ X/Q. For a s determinations, analyt¬ 
ical functions of these power corrections were used to complement the purely 
perturbative expansion, but the more widespread approach to correct for 
non-perturbative effects has been to use Monte Carlo simulations and their 
hadronization models in order to calculate the event shape at the partonic 
and hadronic level. As can be seen in Fig. [3] (right), these non-perturbative 
corrections can be sizeable, especially when approaching the two-jet region of 
the distribution, therefore the fit range has to be chosen carefully. Ultimately 
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the model dependence of such corrections leads to systematic uncertainties 
on the extracted a s values. 


Thrust 



L= — In (1 —T) 



Figure 3: Left: Thrust distribution measured by the ALEPH experiment at 
LEP for centre-of-mass energies between 91.2 and 206 GeV together with 
QCD predictions at NLO plus next-to-leading-logarithmic approximation 
(NLLA). The plotted distributions are scaled by arbitrary factors for pre¬ 
sentation; taken from Ref. |I6]. Right: Comparison of ALEPH data for the 
Thrust distribution to the fitted QCD prediction (NLO+NLLA) obtained at 
parton level (solid line) and corrected for hadronization effects (shaded band). 
The width of the band covers the predictions using different hadronization 
models; taken from Ref. ra- 

A completely different approach to the determination of cc s (M|) consists 
in calculating QCD predictions on the lattice for observables such as hadron 
mass splittings. From the comparison of data to the predictions, obtained as 
a function of the lattice spacing a and extrapolated to a —> 0, first a lattice 
coupling is extracted which is then converted to its perturbative counter-part 
cc s (M|). During the last decades there has been enormous progress in this 
held; indeed, the most precise a s (M|) determinations to date are obtained 
from lattice QCD, though it is fair to say that the community still intensively 
discusses the validity of the very small systematic uncertainties, claimed by 
some of the involved groups. A more detailed discussion of this approach can 
be found in the review by Sachrajda [T8| . 
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4 Brief historical overview 


In the following I make an attempt to sketch some of the relevant develop¬ 
ments that occurred during the last few decades, without any claim of being 
comprehensive and of covering all types of a s studies in a balanced man¬ 
ner. In fact, a particular focus is put on results obtained by experiments at 
CERN. 

The first extensive overview of a s measurements was given by Altarclli 
[19] in 1989. In that review, he summarised measurements based on observ¬ 
ables such as (i) the total hadronic cross section in e + e“ annihilation (at 
that time known at NNLO in pQCD, i.e., up to from data in the range 
7 < Q < 56 GeV; (ii) scaling violations in deep inelastic lcptoproduction with 
structure function data from BCDMS, EMC and CHARM; (iii) quarkonium 
decays, especially ratios of T partial widths (r ggg /r MM , r 7gg /r ggg ) measured 
by the CUSB, CLEO, ARGUS and Crystal Ball collaborations; and (iv) jet 
production, energy-energy correlations and the photon structure function 
from 77 reactions, using e + e“ data collected by the PEP/PETRA experi¬ 
ments. A summary of these measurements is shown in Fig. [4] In an attempt 
to combine all those results, obtained at Q values from a few up to sev¬ 
eral tens of GeV, and using the QCD prediction for the running of a s (Q 2 ), 
he arrived at a prediction for the coupling evaluated at the Z boson mass, 
a s (Q ~ Mz) ~ 0.11 ± 0.01; that is, a determination of the strong coupling 
constant at the 10% level. Interestingly, he concluded with the following 
statement: Establishing that this prediction is experimentally true would be a 
very quantitative and accurate test of QCD, coiiceptually equivalent but more 
reasonable than trying to see the running in a given experiment. It is impres¬ 
sive to note that his prediction turned out to be nicely consistent with the 
current WA value [2] of a s (M|) = 0.1185 ± 0.0006. In addition, we see that 
over the past 25 years the relative uncertainty has been reduced by a factor 
of 18, gauging the enormous progress made during these decades. 

The year 1989 also saw the start of the LEP experiments (ALEPH, DEL¬ 
PHI, L3, OPAL), and the following decade was characterised by great ad¬ 
vances, both experimentally and theoretically, in the field of pQCD in general 
and a s measurements in particular. Extensive overviews can be found, for 
instance, in Refs. |2T, .20] which also discuss the application to earlier JADE 
data of the developments that occurred during the LEP era. 

Event-shape observables were studied in great detail by the LEP exper¬ 
iments, showing first that pQCD at NLO does not provide an adequate ac- 
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Figure 4: Summary of the a s determinations by Altarelli [[jjj] in 1989. 

curacy in order to go well below the 10% level in relative uncertainty on 
a s (M|). At the same time, calculations that resum large logarithms of the 
event-shape variable to all orders in a s appeared and were used to improve the 
NLO prediction for a number of event-shape observables. This also triggered 
the development of a novel jet algorithm, the so-called Durham algorithm 
[22], with a modified metric compared to the previously used JADE algo¬ 
rithm [23]. The modification of the jet metric, which defines the distance in 
phase-space between two particles that might or might not be combined into 
a new pseudo-particle, was motivated by the fact that the pQCD predictions 
for jet rates and the differential 3-jet rate, based on the JADE metric, did not 
show the exponentiation behaviour of the large logarithms 1231, whereas using 
the Durham metric led to exponentiation and ultimately to improved resum¬ 
mation predictions. Note that this Durham algorithm became the standard 
algorithm for jet Ending at LEP, and was at the basis for later develop¬ 
ments of iterative recombination algorithms, nowadays widely used at the 
LHC (cf. Ref. j2]J and references therein). As a consequence of the combined 
NLO+resummed predictions, corrected for non-perturbative hadronization 
effects using phenomenological Monte Carlo (MC) models, the relative un¬ 
certainty of a s was reduced to the ~ 5% level, still dominated by theoretical 
uncertainties due to missing higher orders and estimated from variations of 
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the renormalisation scale. Attempts to replace the MC models by analytical 
power corrections ra of order X/Q did not lead to substantially different 
results. It became clear that only a complete NNLO calculation for jet rates 
and event shapes might lead to a significant reduction of uncertainties. In¬ 
deed, such a calculation [SB:, GZ] appeared after the end of LEP, and its first 
applications [SU to the 3-jet rate and to event shapes, including next-to- 
leading log resummations, resulted in more precise a s measurements at the 
2-3% precision level. 

Most of the aforementioned determinations gave a s (M|) values in a range 
of, very roughly speaking, 0.117 - 0.125. However, as summarized in Ref. [S], 
more recent re-analyses of the Thrust distribution, based on novel develop¬ 
ments in soft-collinear effective held theory, resummation at next-to-next- 
to-next-to-leading logarithmic accuracy, and analytic calculations of non- 
perturbative effects, resulted in values as low as 0.1131, while at the same 
time claiming very small uncertainties at the 2% level. Thus, further work 
will be required to understand this spread of results from jet and event-shape 
observables, which covers a larger range than most of the individually quoted 
uncertainties. 

In terms of inclusive observables, the LEP experiments quoted precise 
a s measurements by using the total hadronic cross section (or equivalently, 
the hadronic width of the Z boson), as well as by analyzing hadronic r 
decays. Contrary to event shapes, NNLO predictions for these observables 
were already available in the nineties, leading to rather small renormalisation 
scale uncertainties. By now, they are even known to N 3 LO accuracy. This 
implies an almost negligible theoretical uncertainty in the case of the hadronic 
Z decay width; for instance, when included in a global fit [3U] to electroweak 
precision data a value of a s (M|) = 0.1193 ± 0.0028 is found, where the 
dominant part of the uncertainty is of statistical nature. 

Naively speaking, similarly precise results might not have been expected 
from the analyses of hadronic r decays, since here the relevant scale is the 
r mass, close to the scale where pQCD is supposed to break down. Thus, 
non-perturbative effects and missing higher order contributions should sig¬ 
nificantly impact the attainable precision. However, it was realised that 
measuring different moments of the r spectral function allows to determine 
a s (M .(?) and to constrain non-perturbative power-suppressed contributions at 
the same time. Several methods were developed to estimate missing higher- 
order terms, beyond NNLO and N 3 LO, such as so-called contour-improved 
perturbative expansions, claiming very small scale uncertainties. It is worth 
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noting [2] that these methods are still matter of intense discussions, in par¬ 
ticular since some of the results obtained by different groups differ by several 
standard deviations in terms of the quoted uncertainties. In an attempt to 
combine all these results and to take into account the observed spread, Ref. 
[2] quotes q; s (M|) = 0.1197±0.0016. Thus, somewhat surprisingly, the study 
of r decays results in one of the most precise a s determinations, basically 
at the level of 1% relative uncertainty. It should be emphasized that the 
already precise measurements, obtained at the scale of the r mass, turn into 
this even more precise result at the Z mass, because of the running of a s as 
discussed at the end of Section [2] 

Many of the developments of the LEP area, in the held of event shapes and 
jet observables, were also applied to HERA data of deep inelastic electron- 
proton scattering (DIS). Although here the pQCD predictions are only known 
up to NLO approximation so far, and the a s extraction from jet final states is 
somewhat more complicated because of the additional implicit a s dependence 
of the PDFs, it is impressive to see that a combination |3Ij based on precise 
HERA data of inclusive jet cross sections in neutral current DIS at high Q 2 
results in a s (M|) = 0.1198±0.0032, which includes a theoretical uncertainty 
of ±0.0026. These HERA measurements also allow covering a large range of 
Q 2 values and thus probing directly the running of a s . 

More inclusive DIS observables, in particular structure functions and their 
scaling violations as discussed in Section [3j have been used in global PDF 
fits based on NNLO pQCD and resulted in smaller relative uncertainties, 
even at the 1% level as quoted by some groups. However, quite similarly 
to the case of event shapes and r observables, also here a spread of a s (M|) 
values (roughly covering a range of 0.113 [32] to 0.117 [33J [31]) is observed 
[2] that is larger than some of the individually quoted uncertainties. Two 
remarks are in place here: (i) these differences are still matter of intense 
discussions among the various groups performing global PDF fits, and (ii) it 
is kind of a tradition that a s determinations from DIS and global PDF fits 
result in smaller values than those obtained from e + e~ annihilations, without 
understanding the origin of this apparent bias. 

Jet observables at hadron colliders, such as the inclusive jet cross section 
as a function of jet transverse momentum or invariant multi-jet masses, jet 
angular correlations or jet rates, are only known to NLO approximation so 
far. Furthermore, important systematic uncertainties due to the jet energy 
scale and choice of PDF set are expected to limit the attainable precision, 
and similarly to the case of DIS, the intrinsic a s dependence of the PDFs 
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has to be carefully taken into account in any a s fit. As discussed in Ref. [2], 
first measurements at the Tevatron and at the LHC gave results consistent 
with the WA value and with relative uncertainties in the range of 4 to 8%. 
However, very important developments have taken place at the LHC recently, 
as e.g. summarised in Ref. |3B], First, in both the ATLAS and CMS experi¬ 
ments the jet energy scale uncertainty is now known at an impressive level of 
about 1-2% for jets in the ~ 100 GeV range [36] [37]. Since jet cross sections 
are steeply falling functions of jet momentum this has an enormous impact 
on the finally attainable precision. Furthermore, ratios of observables, such 
as the ratio of the 3-jet over the 2-jet rate, allow to eliminate this systematic 
uncertainty to a large extend, as shown in Refs. [38, [39]. Finally, NNLO 
calculations for the inclusive jet cross section appear to be around the corner 
m, which will further boost the importance of such measurements. 

In fact, recently the first a s determination j41| at a hadron collider, using 
pQCD at NNLO, has been published. However, here an inclusive quantity, 
namely the top-pair production cross section, has been successfully exploited 
thanks to its strong sensitivity to both a s and the top quark mass. Fixing the 
latter to its WA value allowed the CMS collaboration to extract a s (M|) at an 
impressive relative precision of ~ 3%, also thanks to the remarkable experi¬ 
mental precision (4%) of the top cross section measurement [42] that served 
as input. Because of this recent progress, and because of the large Q 2 range 
covered by the measurements at the LHC, the running of the strong coupling 
constant is now being precisely studied over an unprecedented energy range. 

As mentioned at the end of Section [3] a discussion of a s determinations 
using lattice QCD can be found in a separate review [18] in this volume. For 
completeness it should be stated here that this very complementary approach, 
compared to the measurements described above, results in the world’s most 
precise ay(M|) determinations to date, with some of the analyses quoting rel¬ 
ative uncertainties at the 0.5% level (cf. Ref. [2]). However, the community 
is having intense discussions about the validity of these apparently rather 
optimistic estimates of systematic uncertainties. In any case, the lattice re¬ 
sults dominate the current WA value: not including them in the averaging 
procedure results in cc s (M|) = 0.1183 ± 0.0012 [2], i.e., the uncertainty dou¬ 
bles. 

This historical overview can not be concluded without a brief discussion 
of the general issue of combining a s determinations. As discussed in Ref. [2], 
this is a non-trivial exercise. Since most of the individual measurements are 
dominated by systematic uncertainties, which cannot be expected to follow 
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a normal distribution, and since very often the correlations among these un¬ 
certainties are not very well known, simple averaging methods as applicable 
to measurements with statistical errors only might not be appropriate. In 
1995 Schmelling |43] proposed a method for estimating such unknown corre¬ 
lations, which rescales individual uncertainties according to the assumption 
that the normalized y 2 of the combination should be 1. This method is 
also used for the current WA determination [[2]. Furthermore, there is a cer¬ 
tain arbitrariness in the choice of results included in the average. Finally, 
as mentioned earlier, often a s determinations based on the same observable 
but using different methods give results that differ by a larger amount than 
would be expected from the individually quoted uncertainties, rendering the 
estimate of the combined uncertainty a difficult exercise. 



year 

World average a s (M%) 

Altarelli [19j 

1989 

0.11 ±0.01 

Hinchcliffe fHJ (PDG) 

1992 

0.1134 ±0.0035 

Hinchcliffe [Taj (PDG) 

1995 

0.118 ±0.003 

Schmelling |46] 

1997 

0.118 ±0.003 

Bethke [9] 

2000 

0.1184 ±0.0031 

Bethke [3S] 

2006 

0.1189 ±0.0010 

Bethke, Dissertori, Salam (PDG) [2] 

2013 

0.1185 ±0.0006 



Year 


Figure 5: A selection of world average values for a s (M|) as a function of time; 
the yellow band indicates the current world average value [2] of a s (M £) = 


0.1185 ±0.0006. 
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Throughout these years, several individuals and/or groups have compiled 
the available a s measurements and combined them into a single value. The 
earliest attempt by Altarclli has already been discussed above. During the 
nineties, the reference in terms of a s (M|) was established by the PDG, 
in particular thanks to the PDG review on QCD by Hinchcliffc (see, e.g., 
Refs, m 3Sj)- An independent estimate of the WA value was published by 
Schmelling p|6j in 1997, based on his proposal for handling unknown cor¬ 
relations. Then, during the first decade of this century, Bethke mmm 
provided a number of comprehensive studies, that established the de-facto 
WA value, despite the PDG still publishing an independent combination. 
Since a few years this situation has been resolved, with Bethke now being 
co-author (together with Dissertori and Salam) of the PDG review on QCD 
that also contains the WA determination of a s . Figure [5] displays this, most 
likely incomplete, collection of WA results as a function of time, nicely show¬ 
ing the impressive progress made throughout the last decades. Finally, Fig. [6] 
presents an example [9] of inputs to the averaging procedure and the current 
experimental status of the running of a s , showing excellent agreement with 
the theoretical expectation. 


DIS [pol. strct. fctn.] 
DIS [Bj-SR] 

DIS [GLS-SR] 
x-decays |LEPJ 
xF 3 [v -DIS] 

F 2 [e-, p-DIS] 
jets & shapes [HERA] 
QQ + lattice QCD 
Y decays 



t^e [a had ] .— 

e + e _ [jets & shapes 22 GeV] 
e + e _ [jets & shapes 35 GeV] 

e*-e-[Ohad] 

e + e [jets & shapes 44 GeV] 
e f e _ [jets & shapes 58 GeV] 



pp —> bb X 
PP, PP “> Y X 
o(pp —> jets) 



r(Z°— > had.) [LEP] 
e^e - [scaling. viol.] 
jets & shapes 91.2 GeV [LEP] 



jets & shapes 133 GeV 
jets & shapes 161 GeV 
jets & shapes 172 GeV 
jets & shapes 183 GeV 
jets & shapes 189 GeV 



0.08 0.10 0.12 0.14 

n,(Mz) 



0.1 


Sept. 2013 

x decays (N 3 LO) 

Lattice QCD (NNLO) 

DIS jets (NLO) 

Heavy Quarkonia (NLO) 
e + e“ jets & shapes (res. nnlo) 
Z pole fit (N 3 lo) 
ph -> jets (NLO) 


= QCD a s (M z ) = 0.1185 ± 0.0006 


10 


Q [GeV] 


100 


1000 


Figure 6: Left: List of individual a s (M|) measurements and their comparison 
to the world average from Ref. [9] in 2000; Right: current status of the 
running of a s , as summarised in Ref. J2j. 
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5 Conclusions 


The strong coupling constant is one of the fundamental parameters of the 
standard model of particle physics. In this article I have reviewed the theo¬ 
retical and experimental developments that have led to a precise knowledge 
of this important parameter, representing a cornerstone in our understanding 
of the strong interactions sector of the standard model. 
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